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We report on NMR studies of a honeycomb lattice antiferromagnet Na2Ni2Te06 by ^^Na 

nuclear spin-echo techniques. The ^^Na nuclear spin-lattice relaxation rate l/^^Ti exhibits critical 
divergence near a Neel temperature Tn = 26 K, a narrow critical region, and a critical exponent 
w = 0.34 in 1/^^Ti oc {T/Tjq - 1)“™ for Na2Ni2Te06, and Tn = 18 K for Na2(Nio.5Cuo.5)2Te06. 
Although the uniform magnetic susceptibility of Na2Ni2Te06 exhibits a broad maximum at 35 K 
characteristic of low dimensional spin systems, the NMR results indicate three dimensional critical 
phenomenon around the Neel temperature. 


PACS numbers: 

I. INTRODUCTION 

Na2Ni2Te06 is a quasi-two dimensional honeycomb 
lattice antiferromagnet [ll-[3- The crystal structure 
of Na2Ni2Te06 consists of the stacking of Na and 
(Ni/Te)06 layers (PGs/mcm) ii- The Neel tempera¬ 
ture of Tn ~ 27 K was estimated from measurements of 
specific heat and the derivative of uniform magnetic sus¬ 
ceptibility 0 . The magnetic susceptibility takes a broad 
maximum around 34 K ii- The Weiss temperature of 
0 = - 32 K and the superexchange interaction of J/Zcb = 
- 45 K were estimated from the analysis of Curie-Weiss 
law fit and a high temperature series expansion Q. Al¬ 
though the Ni^+ ion must carry a local moment of 5 = 
1 on the honeycomb lattice, the large effective moment 
/ieff = 3.446/iB could not be explained by spin S = 1 
with a o-factor of g = 2 Q . The g-factor must be larger 
than 2 y , or a Ni^+ ion and the intermediate state might 
be realized because of the tunable valence of Te"’^’*' and 
Te®+ 0. 

Spin frustration effects on a honeycomb lattice have 
renewed our interests, since the discovery of a possible 
spin liquid state in a spin-3/2 antiferromagnet 0. Var¬ 
ious magnetic ground states are competitive with each 
other on the honeycomb lattice 0 . 

In this paper, we report on ^^Na NMR studies of 
Na2Ni2Te06 and Na2(Nio.5Cuo,5)2Te06 polycrystalline 
samples. Na2(Nio,5Cuo.5)2Te06 still beloMS to the same 
space group PQ^/mcm as Na2Ni2Te06 @, 0. For the 
Cu substitution, we expected a possible enhancement of 
quantum effects from S' = 1 to 1/2. Since the solubility 
limit in the honwcomb lattice Na2(Nii_2;Cua;)2Te06 is 
about a; = 0.6 0 , we selected the half Cu-substituted 
sample being away from the phase boundary. We ob¬ 
served three dimensional critical phenomenon in the 
^^Na nuclear spin-lattice relaxation rate 1/^^Ti near 
Tn = 26 K for Na2Ni2Te06 and Tn = 18 K for 
Na2(Nio.5Cuo,5)2Te06. The broad maximum of uni¬ 
form magnetic susceptibility is not the onset of mag¬ 
netic long range ordering. In the antiferromagnetic state 


of Na2Ni2Te06, we observed 1/^^TiocT^ which indicates 
conventional spin-wave scattering. 

II. EXPERIMENTS 

Powder samples of Na2Ni2Te06 have been synthesized 
by a conventional solid-state reaction method. Appropri¬ 
ate amounts of NiO, TeOe and Na2C03 were mixed, pal¬ 
letized and fired a few times at 800 - 860°C and finally at 
900°C for 24 hours in air. The products were confirmed 
to be in a single phase from measurements of powder X- 
ray diffraction patterns. Magnetic susceptibility x at 1.0 
T was measured by a superconducting quantum interfer¬ 
ence device (SQUID) magnetometer. Powder samples of 
Na2(Nio.5Cuo.^2Te06 were previously synthesized and 
characterized 0. 

A phase-coherent-type pulsed spectrometer was uti¬ 
lized to perform the ^^Na NMR (nuclear spin / = 3/2) 
experiments in an external magnetic field of 7.4847 T. 
The NMR frequency spectra were obtained from Fourier 
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FIG. 1: Uniform magnetic susceptibility x of Na2Ni2Te06 
and Na2(Nio.sCuo,5)2Te06. Solid curves are the results from 
least squares fits by a Curie-Weiss law. 
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transformation of the ^^Na nuclear spin-echoes. The 
^^Na nuclear spin-lattice relaxation curves = 1 — 

Eit)/E{oo) (recovery curves) were obtained by an in¬ 
version recovery technique as a function of time t after 
an inversion pulse, where the nuclear spin-echoes E(t), 
E{oo)[= i?(10Ti)] and t were recorded. 


III. EXPERIMENTAL RESULTS AND 
DISCUSSIONS 


least-squares fits by a theoretical multi-exponential func¬ 
tion for a central transition line (/^ = 1/2 0 -1/2) 

23p(i) = p(0){0.1e“‘/”^i -b (1) 

where p{0) and a ^^Na nuclear spin-lattice relaxation 
time are fit parameters. The theoretical function of 
eq. m well reproduces the experimental recovery data. 
Thus, the assignment of the exciting spectrum to the 
central transition line is justified a posteriori^ too. 


A. Uniform magnetic susceptibility 

Figure [T] shows uniform magnetic susceptibility x 
of Na2Ni2Te06 and Na2(Nio.5Cuo,5)2Te06. The solid 
curves are the results from least squares fits by a Curie- 
Weiss law. We estimated the Weiss temperature 9 = - 27 
K and the effective moment peS = 3.4/iB for Na2Ni2Te06, 
which agree with the previous report Q, and 0 = - 35 
K and peS = 2.5^b for Na2(Nio.5Cuo,5)2Te06. If the g- 
factor is 5 = 2, then S' = 1 and S = 1/2 lead to pes = 
2.83/iB and 1.73/iB, respectively, x deviates below about 
100 K from the Curie-Weiss law and takes a broad max¬ 
imum at 35 K in Na2Ni2Te06. x drops below about 20 
K in Na2(Nio,5Cuo.5)2Te06. 


B. NMR spectrum and recovery curves 


(a) (b) 



FIG. 2: (a) Fourier-transformed NMR spectrum at 

84.670 MHz and at 300 K. (b) ^^Na nuclear spin-lattice re¬ 
laxation curves ^^p(t) at a central frequency. Solid curves are 
the results from least squares fits by eq. (HJ. 


Figure[2](a) shows the Fourier-transformed spectrum of 
23Na spin-echoes at a Larmor frequency of 84.670 MHz 
and at 300 K. The central transition line of Iz = 1/2 tA 
-1/2 is affected by a nuclear quadrupole interaction Q. 
The linewidth is about 150 kHz. The precise value of the 
Knight shift could not be determined within the present 
studies. 

Figure[2Kb) shows the recovery curves ^^p(t) with vary¬ 
ing temperature. The solid curves are the results from 


C. Na2Ni2Te06 


(a) 
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FIG. 3: (a) l/^^Ti and uniform magnetic susceptibility y 
against temperature. shows a critical divergence near 

Tn = 26 26.5 K and levels off above about 100 K. The 

broken line indicates l/^^Tioo = 88 s“^. (b) l/^^Ti and y 
against temperature in enlarged scales. Solid curves are visual 
guides. 


Figures IS a) and (b) show 1/^^Ti and uniform mag¬ 
netic susceptibility y against temperature. 1/^^Ti takes 
l/^^Tioo = 88 s“^ above about 100 K and shows a di¬ 
vergence at 26 ^ 26.5 K which can be assigned to the 
Neel temperature Tn. Thus, the broad maximum of the 
magnetic susceptibility y at 35 K is not due to the an¬ 
tiferromagnetic long range ordering but due to a low di¬ 
mensional short range correlation developing on the hon¬ 
eycomb lattice antiferromagnets Q . The result is consis¬ 
tent with the specific heat measurements Q. 
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FIG. 4: (a) 1/^^Ti against temperature. The solid curve is 
the result from a least squares fit by eq. The Neel tem¬ 
perature and the critical exponent were estimated to be Tn 
= 26.24 K and w — 0.34, respectively, (b) Log-log plots of 
normalized (l/^®Ti)/(l/^®Tioo) against reduced temperature 
|T —TnI/Tn. Closed and open triangles are above 

and below Tn , respectively. The solid line indicates the result 
from a least squares fit by eq. ©. 

FigureSlJa) shows 1/^^Ti against temperature and the 
result (the solid curve) from a least-squares fit by 


23Ti 232^100 iT/TN-ir’ 

where a constant (7, a Neel temperature Tn, and a critical 
exponent w are fit parameters. The fitting results were 
Tn = 26.24 K and w = 0.34. 

A mean field theory for a three dimensional isotropic 
Heisenberg antiferromagnet gives tc = 1/2 Q. A dynamic 
scaling theory gives w = 1/3 for a three dimensional 
isotropic Heisenberg model [lOl| and w = 2/3 for a three 
dimensional uniaxial anisotropic Heisenberg model [HI- 
The exponent of rc = 0.34 indicates that Na2Ni2Te06 
in the critical region is described by a three dimensional 
dynamical spin susceptibility. In passing, CuO exhibits 
a similar w = 0.33, a broad maximum in x at 540 K, and 
Tn = 230 K [Si- 

Figure SKb) shows log-log plots of normalized 
(1/23 Ti)/(1/23Tioo) against reduced temperature 


FIG. 5: (a) l/^^Ti against temperature for Na2Ni2Te06 

and Na2(Nio.5Cuo.5)2Te06. The broken lines indicate 
1/^^Tioo = 88 and 57 s“^. (b) Log-log plots of 

normalized against reduced tempera¬ 

ture IT-TnI/Tn for Na2Ni2Te06 (Tn = 26.24 K) and 
Na2(Nio.5Cuo.5)2Te06 (Tn = 18 K). The solid line is eq. Q 
with the critical exponent of ui = 0.34. 

|T — Tn|/Tn. The solid line indicates the result from a 
least squares fit by eq. @. 

The onset of increase in the NMR relaxation rate near 
Tn empirically categorizes critical regions. The region 
of |T — Tn|/Tn < 10 has been assigned to the renor¬ 
malized classical regime with a divergent magnetic cor¬ 
relation length toward T = 0 K [S- The region of 
|T — Tn|/Tn < l.O has been assigned to the three di¬ 
mensional critical regime with a divergent magnetic cor¬ 
relation length toward Tn. Thus, the narrow critical 
region of |T —Tn|/Tn < I also empirically categorizes 
Na 2 Ni 2 Te 06 to the three dimensional critical regime. 

At high temperatures of T ^ J, the spin system is in 
the exchange narrowing limit. Then, l/^^Ti is expressed 

by 

1 /:y^*5'(5'-b 1) 2:n( 7nA) 

_ = ,/2¥^-, (3) 

^T = |s(S + l)i(2)t (4) 
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FIG. 6: Log-log plots of against temperature for 

Na2Ni2Te06 and Na2(Nio.5Cuo.5)2Te06. A broken line in¬ 
dicates a function of eq. ©• Solid curves are visual guides. 


where ^^ 7 n,/ 27 r = 11.262 MHz/T is the nuclear gy- 
romagnetic ratio, A is a hyperfine coupling constant, and 
ujex is an exchange frequency [l^. is the number of 
Ni ions nearby a ^^Na nuclear, z is the number of the 
nearest neighbor Ni ions. Assuming J = 45 K, Q S' = 
1, and z = 3, we obtained Wex = 12 x 10^^ s“^. From 
eq. dni) with 1/^^Tioo = 88 s“^, we derived the hyperfine 
coupling constant A = 2.0 kOe/uu, which is nearly the 
same as that of Na 3 Cu 2 Sb 06 [iq . 


D. Na2(Nio.5Cuo.5)2Te06 

Figure [SJa) shows 1/^^Ti against temperature for 
Na2Ni2Te06 and Na2(Nio.5Cuo.5)2Te06. For the half 
substitution of Cu for Ni, 1/^^Tioo and Tn decrease to 
57 s“^ and 18 K, respectively. Extrapolating linearly 
Tn with ATn = - 8 K per half Cu to full Cu substitu¬ 
tion, one may infer Tn = 10 K of a hypothetical spin- 
1/2 honeycomb lattice Na 2 Cu 2 Te 06 ,” although actual 
Na 2 Cu 2 Te 06 is known to be monoclinic and an alternat¬ 
ing spin chain system Bill- 

Figure [Sl[b) shows log-log plots of normalized 
(l/^^Ti)/(l/^^T'ioo) against reduced temperature 
|T-Tn|/Tn for Na2Ni2Te06 (Tn = 26.24 K) and 
Na2(Nio,5Cuo.5)2Te06 (In = 18 K). The solid line 
indicates eq. @ with the critical exponent oi w = 


0.34. The critical region of Na2(Nio.5Cuo.5)2Te06 is still 
narrow as the same as that of Na2Ni2Te06. Simply, Tn 
decreases. No dimensional crossover is observed. 


E. Below Tn 


Figure [6] shows log-log plots of 1/^^Ti against temper¬ 
ature for Na2Ni2Te06 and Na2(Nio.5Cuo.5)2Te06. With 
cooling down below Tn, 1/^^Ti rapidly decreases. The 
broken line indicates a T^ function as a visual guide. In 
conventional antiferromagnetic states, the nuclear spin 
transitions are caused by Raman scattering and three 
magnon scattering [l^ . Then, 1 /Ti is expressed by 



( 5 ) 


in a temperature range of Tn > T ^ Tae, where 
Tae corresponds to an energy gap in the spin wave 
spectrum [^. The energy gap is due to a crystalline 
anisotropy field. The rapid drop of 1/^^Ti below Tn re¬ 
sults from the suppression of low energy excitations by 
the energy gap. Below Tae, an activation-type temper¬ 
ature dependence should be observed in 1/Ti. Since no 
activation behavior was observed down to 5 K, one may 
estimate Tae < 5 K. 


IV. CONCLUSIONS 

In conclusion, we found three dimensional critical phe¬ 
nomenon near Tn = 26 K for Na2Ni2Te06 and Tn = 18 
K for Na2(Nio.5Cuo.5)2Te06 from measurements of the 
^^Na nuclear spin-lattice relaxation rate 1/^^Ti. We have 
analyzed the NMR results by Ni^+ with S' = 1 and ob¬ 
tained sound values of parameters for Na2Ni2Te06. We 
attribute the deviation from the Curie-Weiss law and the 
broad maximum of uniform magnetic susceptibility to 
two dimensional spin-spin correlation on a honeycomb 
lattice. 

The author thanks M. Isobe (Max Planck Institute) 
for X-ray diffraction measurements, K. Morimoto, C. 
Michioka, K. Yoshimura (Kyoto University) for sample 
preparation and characterization at an early stage. 
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